Tri-critical behavior in rupture induced by disorder 
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Abstract 

We discover a qualitatively new behavior for systems where the load transfer 
has limiting stress amplification as in real fiber composites. We find that the 
disorder is a relevant field leading to tri-criticality, separating a first-order 
regime where rupture occurs without significant precursors from a second- 
order regime where the macroscopic elastic coefficient exhibit power law be- 
havior. Our results are based on analytical analysis of fiber bundle models 
and numerical simulations of a two-dimensional tensorial spring-block system 

in which stick-slip motion and fracture compete. 
05.50.+q, 46.30.Nz, 81.40.Np 
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There is growing evidence that rupture in random media can be viewed as a kind of 
critical phenomenon as a result of the interplay between disorder and fracture mechan- 
ics, with proposed applications in particular to fiber composites |3],|4j] and earthquakes ||. 
Notwithstanding its importance, we do not have a comprehensive understanding of rupture 
phenomena but only a partial classification |TJ,|2|||]. /,From a theoretical point of view, rup- 
ture is controlled in principle by the infinite moment (a q )\ q -+ 00 of the stress field and the 
difficulties emerge from the non-commutation of the two limits (q — > oo, A — > 0), where A is 
the amount of disorder (see below for a precise definition). In intuitive wording, the largest 
stress in the system is very sensitive to the amount and type of disorder. Disorder is known 
to induce stress field distributions with fat tails 0. Consider for instance a log- normal 
distribution with standard deviation A and mean o"o, then (a q ) l l q = <7 e A2q / 2 , which shows 
that the limit A — > is singular for rupture (q — > oo). This non-commutativity of limits 
is at the crux of some of the major outstanding problems in physics || such as turbulence 
(viscosity — > ; time — > oo) and quantum chaos (h — > ; time — > oo). Here, we show that 
the amount of disorder A plays the role of a relevant field which makes systems with limited 
stress amplification exhibit a tri-critical transition as the disorder increases, from a Griffith- 
type abrupt rupture (first-order) regime to a progressive damage (critical) regime. This is 
reminiscent of the critical behavior induced by quenched disorder in magnetic systems ||. 

We first document this behavior in a simple mean-field model of rupture, known as the 
democratic fiber bundle model fllPf . It consists of N parallel fibers with identical spring 
constants and identically independent random failure thresholds Xj distributed according 
to the cumulative probability distribution P{Xj < x) = P(x). A total force F is applied to 
the system and is shared democratically among the N fibers. When the force on one fiber 
reaches its threshold, the fiber ruptures and the stress is redistributed to all remaining fibers. 
This transfer might induce secondary failures which in turn induce tertiary ruptures and so 
on. One is interested in the stress-strain characteristic as the applied force is increased, the 
properties of the rupture point and the precursory events prior to the complete breakdown. 
The solution of this problem is found by noticing that the total bundle will not break under 
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a load F if there are n fibers in the bundle each of which can withstand x n = F/n. x n and 
n are related, for large N, by n = N[l — P(x n )] leading to F(x n ) = Nx n [l — P(x n )]. The 
number k of fibers which have failed under the force F is then k = N — n = NP(x n ). Now, 
for a broad class of distribution P(x) extending down to 0, the function x[l —P(x)] presents 
a maximum at < x* < oo, the solution of dx[l — P{x)\/dx\ x=x * = 0. As the behavior of 
F{x n ) close to x* is quadratic F(x n ) ~ F* — c(x* — x n ) 2 where c is a constant, this implies 
that the rate dk/dF of fiber failure diverges as 

(F*-F)~ 1 /2 ) w h er e F* = x*[l-P{x*)}, thus 

qualifying a critical mean field behavior. 

However, if P(x) is such that dx[l — P(x)]/dx = has no solution, the behavior will 
be completely different, with a sequence of a few fibers maybe breaking as the load is 
applied followed by an abrupt global failure. Correspondingly, the stress-strain characteristic 
exhibits a discontinuity in its slope at the point of rupture. This qualifies a first-order 
behavior. Notice that this is similar to the Ehrenfest's classification of the order of phase 
transitions, where the free energy is here replaced by the elastic energy. 

Let us take for instance P(x) = for < x < xi, P(x) = {x — X\)/ A for xi < x < x± + A 
and P(x) = 1 for x > X\ + A, corresponding to the strengths Xj uniformly distributed 
between x± > Oandxi+A. Then, dx[l— P{x)\/dx = (x\+A—2x)/A, which has a root in the 
interval x\ < x < Xi+A if and only if A > x±. In this case, we recover the previous mean field 
critical behavior. However, for weak disorder A < xi, not a single fiber breaks down until the 
force reaches iVxi at which value the system of N fibers breaks suddenly. This is an extreme 
illustration of a "first-order" behavior. The particular value A = x\, F = Nx± thus plays 
the role of a tri-critical point in analogy with thermal phase transitions ||11|| . This behavior 



holds for a large class of distributions P(x): the condition that dx[l — P(x)]/dx has no root 
is equivalent to the condition that the equation dlog[l — P(x)]/dx — — 1/x has no solutions 
for any < x < oo. This equation defines two domains: (1) dlog[l — P(x)]/dx < —1/x for 
all x > 0: this can occur in particular if 1 — P(x) decays to zero faster than 1/x for large 
x, with the additional constraint that there exists a minimum strength x\ strictly positive. 
Notice that the distributions which extend down to zero are in this sense always in the 



"large" disorder regime. (2) dlog[l — P{x)\/dx > — 1/x for all x > 0: this corresponds to 
distributions which decays slower than 1/x. Take for instance 1 — P(x) = (1 +x)~ a . Then, 
dlog[l — P(x)]/dx = —a/(l + x) which remains strictly larger than —1/x if a < 1. 

Having established the existence of the tri-critical behavior in this mean field model, 
let us now turn our attention to a more realistic two-dimensional (2D) spring-block model 
of surface fracture in which the stress can be released by spring breaks and block slips. 
We consider the experimental situation where a balloon covered with paint or dry resin is 
progressively inflated . An industrial application is a metallic tank with carbon or kevlar 
fibers impregnated in a resin matrix wrapped up around it which is slowly pressurized [[|. 
As a consequence, it elastically deforms, transferring tensile stress to the overlayer. Slipping 
(called fiber-metal delamination) and cracking can thus occur in the overlayer. We model 
this process by a 2D array of blocks which represents the overlayer on a coarse grained 
scale in contact with a surface with solid friction contact. The solid friction will limit stress 
amplification. Each block is interconnected to its nearest neighbors via springs of unstretched 
lengths Iq and spring constants K. The position of the blocks in the x- and ^-directions 
are given by (a • % + Xjj, a • j + y^j) where 1 < i, j < L, form a square lattice with lattice 
constants a, and where Xij, yij fulfill x^^y^ <C a, so that Hooke's law applies. In |TS| it 



was shown that the x component of the force on a block to first order in the displacements 
takes the form: 

Fij — ~K{{bi + \j + bi_ij)xij — bi + ijx i+ ij — 

bi-i,jXi-i,j + s [(bi,j+i + bij-ijXij — bij-iXij-i 

-b id+1 x iJ+ i] - as(b i+1J - fei-ij)} (!) 

and, by symmetry, Ffj follows by switching x y and i <-> j. s = (a — lo)/a > is the 
strain of the network without fluctuations (xij,yij = 0) , and &j±i,j±i = 1,0, respectively, 
depending on whether a spring connects the blocks (i, j) — (i ± 1, j ± 1) or not. Likewise the 
stress B in a spring is given by: 

B(i,j)-{i±l,j) = K[{ x i,j — x i±l,j — s ) 2 + 
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s 2 (y i , j -y i ±i, j ) 2 ] 1/2 (2) 

Initially Xij and yij are chosen uniformly from the interval [— A, +A], thus A quantifies 
the amount of disorder which is on the initial displacements corresponding to an effective 
initial disorder in the thresholds. Periodic boundary conditions are used in both the x- and 
^-direction. 

The coupling of the overlayer to the substrate has two effects when the substrate expands: 
(1) Tensile stress is transferred to the overlayer. This is taken into account by imposing an 
increase in the average distance a between the blocks so as to reflect the inflation of the 
balloon. As a definition of the time, t, we let a(t) increase linearly with t. (2) The increasing 
tensile stress in turn gives rise to stick-and-slip motion or/and cracking. A block is assumed 
to stick until the total force applied on it exceeds a threshold F s , where after it slips to the 
zero-force position, corresponding to local mechanical equilibrium in the absence of friction. 
This thereby releases stresses on its neighbor blocks. A spring breaks irreversibly once the 



stress B exceeds a threshold F r = kF„ 14 



We define a time dependent apparent macroscopic stress on the system, a app {t), from 
the relation a app (t) = E(t)/e(t), where E(t) is the total elastic energy stored in the springs 
of the system at time t, and e(t) = a(t)/a(t = 0) is the macroscopic strain. We calculate an 
effective Young modulus, given by 

Y app (t) = da app (t)/de. (3) 

Y app (t) can be expressed as [da app (t)/da(t)](da/de), where a is proportional to the first 
invariant (the trace) of the real stress field in the system, da/de is the corresponding 
elastic modulus expected to exhibit a power law behavior if criticality is present, while 
da app (t)/da{t) goes to a constant. Therefore, the measurement of Y app (t) gives us direct 
access, if present, to the critical behavior of the Young modulus of the system. As the strain 
e(t) is increased, block slips and spring failures occur up to a point where the system is 
completely ruptured and the stress necessary to impose a constant small strain rate starts 
to decrease from a maximum. At this point, there is at least one large crack spanning the 



whole system. If global rupture occurs abruptly (first-order case), cr app (t) must exhibit a 
sharp maximum and Y app (t) remains finite. If, on the other hand, the rupture is critical, 
Cappit) will exhibit a progressive rounding with a smooth maximum, while Y app (t) vanishes 
as a power law Y app oc [(e c — e)/e c ] 7 on the approach to rupture at e c . 

In Fig. 1 are shown the stress-strain curves from one single realization for each different 
system size, with A = 0.75A C , where A c is the maximal amplitude such that at t — 
B(i,j)-(i±i,j±i) < F c an d < F s for ah (hj)- A c has been determined numerically for 
each run. We set a(t = 0) = lo = 1 and Fs/K = 1 throughout this paper. Observe that the 
maximum stress a system can sustain is an increasing function of re and the range of e- values 
over which fracturing takes place decreases as re increases. For A = 0.75A C and re < 2.9, the 
stress-strain curve presents a smooth maximum indicating a critical rupture. This is con- 
firmed in Fig. 2, showing the vanishing of the apparent Young modulus Y app as e — > e c from 
below. Each curve is obtained by averaging over N = 1000 — 5000 independent configura- 
tions with system size L = 30. e c has been estimated from the condition d{a app ) /de\ e=€c = 0, 
where (• • •) stands for an ensemble average. Larger lattice sizes (L = 50 — 400) with the 
same value of LN were used, but for a given fixed LN we found the smallest lattice sizes L 
(= 30) give the best statistics, which we attribute to the lack of self averaging. For small re, 
the exponent 7 approaches a value slightly larger than 1, while it decreases continuously to 
zero as re increases, as shown in the inset. It seems to vanish around re = 2.9, signaling the 
transformation of the critical regime to an abrupt "first-order" behavior. Keeping re fixed 
and varying A/A c , we find that 7 stays constant, but the size of the critical region increases 
with the magnitude of A/A c . It shrinks to zero at a threshold value function of re which 
is shown in Fig. 3. This function gives the boundary in the (A/A c ; re) plane between the 
critical and first-order regime. As announced, for fixed re < 2.9, increasing the disorder A 
allows the system to go from a first-order to a critical regime. The fact that the disorder 
is so relevant as to create the analog of a tri-critical behavior can be tracked back to the 
existence of solid friction on the blocks which ensures that the elastic forces in the springs 
are carried over a bounded distance (equal to the size of a slipping "avalanche" ) during the 
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stress transfer induced by block motions. 

When k is large, the system responds initially to an expansion by the release of stress 
uniquely by block slips. The block slips give rise to a stress rearrangement, and a spatial 



coarsening phenomenon of the stress field B takes place |T3| . The enhanced correlations in 
the .B-field result in turn in a coherence when fracturing sets in, amounting to smoothing 
out the disorder, thus allowing for a large crack to develop in an abrupt way. For sufficiently 
large k, the system breaks into two parts. This is the regime of first-order behavior. Notice 
that increasing k in the (A/A c ; k) phase diagram corresponds to decreasing the disorder and 
changing at the same time the distribution of disorder so that it becomes more correlated. 
This is therefore a more complicated route than just changing the width of the threshold 
distribution as in the previous fiber bundle model. 

The value of the Young critical exponent 7 for small k can be predicted from percolation 
theory. Indeed, consider the limit k — > 0, for which the blocks are stuck to the substrate and 
cannot move. Only the springs can fail and they do so in a completely uncorrelated way, 
controlled by the initial random configuration of the blocks. We thus get an uncorrelated 
random dilution, ending at the percolation threshold where a macroscopic crack spans the 
system. In the presence of internal strain, it was shown that the elastic constant decreases 
to zero when the dilution increases with an exponent given by the scalar elasticity problem 
||16||, equal to the conductance exponent of percolation. Extensive numerical simulations 



give the value 1.300 for this exponent ||17|| , which is in agreement with the extrapolation of 
our results for k — > 0. 

It is important to understand that these properties belong to systems with load transfer 
mechanisms limiting stress amplification at crack tips. If no coupling or delay mechanism 
exist to reguralize the divergence induced by elasticity at the crack tips (with a stress 
diverging as l/y/r next to a crack tip), the first-order behavior is only observed for zero 
disorder described by the single-crack Griffith criterion and any amount of disorder is relevant 
to produce a critical behavior JI|Q]. However, even in this case, the amount of disorder 
remains of utmost importance as it controls the size of the critical region, and therefore its 



observability |L8 



The existence of different regimes for rupture, depending on the limiting stress amplifica- 
tion and on disorder, opens the road to important potential applications for failure prediction 
purposes such as in the time-to-failure approach fll9~| . We suggest that the often observed 
power law distribution of acoustic emission bursts of many materials upon stressing, offers 
an additional evidence of the critical nature of the damage and cracking of heterogeneous 
materials . Our results provide the foundation for understanding why some systems 
exhibit clearer precursors before rupture than others in which they may even be absent in 
certain cases and for quantifying the expected amount and style of precursory activity as a 
function of heterogeneity and range of interaction. 
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FIGURES 

FIG. 1. Stress-strain curves for different values of k and for different system sizes L = 50 
(dotted line), 100 (thin bold line), and 200 (bold line). The fracturing is stopped at different e(t) 
for different L in order to distinguish between the curves. One L = 400 simulation has been done 
for k = 4 (fat bold line). 

FIG. 2. Macroscopic Young modulus vs reduced macroscopic strain for different values of 
k = 0.5(o), 0.75(+), 0.875(D), 1.4(x) and 2.0(A). A/A c = 0.75. The inset shows the exponent 7 
as a function of k. k « 0.5 is the smallest value for which the system initially has no bonds that 
exceed the threshold. 



FIG. 2. Inset to Fig. 2 

FIG. 3. Phase diagram for criticality of the fracturing. A/A c — > can not been studied 
within the spring-block model since ambiguity in the updating rules for stress release, affects the 
fracturing. 
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